In [2]:

    

using DataFrames
df = readtable("../../results/shadowprice.csv")


    

    
Out[2]:




x1
x2
fips
1
0.0
0.0
1001
2
0.0
0.0
1003
3
0.0
0.0
1005
4
0.0
0.0
1007
5
0.0
0.0
1009
6
-3.78945111378095
-0.8692009518975058
1011
7
0.0
0.0
1013
8
0.0
0.0
1015
9
-3.7345764162987303
-0.8692009518975058
1017
10
0.0
0.0
1019
11
0.0
0.0
1021
12
0.0
0.0
1023
13
0.0
0.0
1025
14
0.0
0.0
1027
15
0.0
0.0
1029
16
0.0
0.0
1031
17
0.0
0.0
1033
18
0.0
0.0
1035
19
0.0
0.0
1037
20
0.0
0.0
1039
21
0.0
0.0
1041
22
0.0
0.0
1043
23
0.0
0.0
1045
24
0.0
0.0
1047
25
0.0
0.0
1049
26
0.0
0.0
1051
27
0.0
0.0
1053
28
0.0
0.0
1055
29
0.0
0.0
1057
30
0.0
0.0
1059
&vellip
&vellip
&vellip
&vellip

In [3]:

    

include("../../src/lib/graphing.jl")
include("../../src/lib/datastore.jl")


    

    




INFO: Precompiling module RCall...
ERROR: LoadError: LoadError: Could not find R installation. Try setting "R_HOME" environmental variable.
 in locate_rhome at /Users/jrising/.julia/v0.4/RCall/src/setup.jl:83
 in include at /Applications/Julia-0.4.0.app/Contents/Resources/julia/lib/julia/sys.dylib
 in include_from_node1 at /Applications/Julia-0.4.0.app/Contents/Resources/julia/lib/julia/sys.dylib
 in include at /Applications/Julia-0.4.0.app/Contents/Resources/julia/lib/julia/sys.dylib
 in include_from_node1 at /Applications/Julia-0.4.0.app/Contents/Resources/julia/lib/julia/sys.dylib
 [inlined code] from none:2
 in anonymous at no file:0
 in process_options at /Applications/Julia-0.4.0.app/Contents/Resources/julia/lib/julia/sys.dylib
 in _start at /Applications/Julia-0.4.0.app/Contents/Resources/julia/lib/julia/sys.dylib
while loading /Users/jrising/.julia/v0.4/RCall/src/setup.jl, in expression starting on line 86
while loading /Users/jrising/.julia/v0.4/RCall/src/RCall.jl, in expression starting on line 21






    




LoadError: LoadError: Failed to precompile RCall to /Users/jrising/.julia/lib/v0.4/RCall.ji
while loading /Users/jrising/research/water/awash/src/lib/graphing-rcall.jl, in expression starting on line 1
while loading /Users/jrising/research/water/awash/src/lib/graphing.jl, in expression starting on line 2

 in error at /Applications/Julia-0.4.0.app/Contents/Resources/julia/lib/julia/sys.dylib
 in compilecache at loading.jl:383
 in require at /Applications/Julia-0.4.0.app/Contents/Resources/julia/lib/julia/sys.dylib
 in include at /Applications/Julia-0.4.0.app/Contents/Resources/julia/lib/julia/sys.dylib
 in include_from_node1 at /Applications/Julia-0.4.0.app/Contents/Resources/julia/lib/julia/sys.dylib
 in include at /Applications/Julia-0.4.0.app/Contents/Resources/julia/lib/julia/sys.dylib
 in include_from_node1 at /Applications/Julia-0.4.0.app/Contents/Resources/julia/lib/julia/sys.dylib

In [11]:

    

include("../../src/lib/readconfig.jl")
config = readconfig("../configs/standard-1year.yml")
suffix = getsuffix()


    

    
Out[11]:





""

In [12]:

    

usmap(DataFrame(fips=df[:fips], value=df[:x1]))


    

    




WARNING: RCall.jl: Warning: package ‘ggplot2’ was built under R version 3.2.5
Use suppressPackageStartupMessages() to eliminate package startup
messages.
WARNING: RCall.jl: 
-----------------------------------------------------------
PBS Mapping 2.69.76 -- Copyright (C) 2003-2017 Fisheries and Oceans Canada

PBS Mapping comes with ABSOLUTELY NO WARRANTY;
for details see the file COPYING.
This is free software, and you are welcome to redistribute
it under certain conditions, as outlined in the above file.

A complete user guide 'PBSmapping-UG.pdf' is located at 
/Library/Frameworks/R.framework/Versions/3.2/Resources/library/PBSmapping/doc/PBSmapping-UG.pdf

Packaged on 2015-04-23
Pacific Biological Station, Nanaimo

All available PBS packages can be found at
http://code.google.com/p/pbs-software/

To see demos, type '.PBSfigs()'.
-----------------------------------------------------------


WARNING: RCall.jl: Loading required package: maptools
Loading required package: sp
Warning: package ‘sp’ was built under R version 3.2.5
Checking rgeos availability: TRUE
Loading required package: foreign






    













    
Out[12]:





RCall.RObject{RCall.VecSxp}

In [13]:

    

df = readtable("../results/shadowprice-1000.csv")


    

    
Out[13]:




x1
x2
fips
1
-100.0
-100.0
1001
2
-100.0
-99.99999999999999
1003
3
-100.0
-100.0
1005
4
-100.0
-99.99999999999999
1007
5
-100.00000000000001
-100.0
1009
6
-100.0
-100.0
1011
7
-100.0
-100.0
1013
8
-100.0
-100.0
1015
9
-100.0
-100.0
1017
10
-100.0
-100.0
1019
11
-100.0
-99.99999999999999
1021
12
-100.0
-100.0
1023
13
-100.0
-100.0
1025
14
-100.0
-100.0
1027
15
-100.00000000000001
-100.00000000000001
1029
16
-100.0
-100.0
1031
17
-100.0
-100.00000000000001
1033
18
-100.00000000000001
-100.0
1035
19
-100.0
-100.0
1037
20
-100.0
-99.99999999999999
1039
21
-100.0
-100.0
1041
22
-100.00000000000001
-99.99999999999999
1043
23
-100.0
-100.00000000000001
1045
24
-100.0
-100.0
1047
25
-100.0
-100.00000000000001
1049
26
-100.0
-100.00000000000001
1051
27
-100.0
-100.0
1053
28
-100.00000000000001
-100.00000000000001
1055
29
-100.0
-100.0
1057
30
-100.0
-100.00000000000001
1059
&vellip
&vellip
&vellip
&vellip

In [14]:

    

usmap(DataFrame(fips=df[:fips], value=df[:x1]))


    

    













    
Out[14]:





RCall.RObject{RCall.VecSxp}

In [15]:

    

df = readtable("../results/shadowprice-1e6.csv")
usmap(DataFrame(fips=df[:fips], value=df[:x1]))


    

    













    
Out[15]:





RCall.RObject{RCall.VecSxp}

In [ ]: